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ABSTRACT 

The  objective  of  this  paper  is  to  present  a  new  computational  methodology  for  predicting  the  durability  of 
structures  made  of  composite  materials  based  on  epoxy  matrix  with  long  carbon  fibres. 

To  analyse  the  behaviour  of  composite  materials  an  Enhanced  Serial-Parallel  constitutive  model  (ESP  model) 
is  developed  assuming  components  behave  as  parallel  materials  in  the  fibres  alignment  direction  and  as  serial 
materials  in  orthogonal  directions.  It  allows  equal  component  strains  in  the  fibre  direction  and  equal  stresses 
in  the  transverse  directions.  The  ESP  model  brings  answers  on  the  non-linear  behaviour  of  composites, 
where  Classical  Micro-mechanics  Formulas  are  restricted  to  their  linear  elastic  part.  Constitutive  tensors  of 
the  composite  materials  are  obtained  from  mechanical  properties  and  internal  variables  of  their  components 
within  a  continuum  framework.  Anisotropy  and  different  constitutive  models  (plasticity,  damage  and  fatigue) 
for  each  phase  are  considered. 

The  fatigue  prediction  model  is  based  in  a  continuum  mechanic  stress  life  approach  for  each  of  the  component 
materials.  S-N  curves  are  proposed  for  each  phase.  A  cumulative  fatigue  damage  index  is  used  to  update  the 
mechanical  properties  of  the  components  and  with  these  properties  update  layers  properties  and  compute 
residual  stiffness  of  the  laminate. 

Constitutive  models  for  the  composite  materials  have  been  implemented  into  a  general  purpose  finite  element 
code  (COMET)  and  a  user  friendly  interface  for  composite  material  data  input  in  a  pre-post  processing 
module  ( GiD)  has  been  developed.  The  methodology  is  validated  using  experimental  data  from  tests  on  CFRR 
composite  material  samples. 

1.0  INTRODUCTION 

Fatigue  is  defined  as  "the  process  of  permanent,  progressive  and  localized  structural  change  which  occurs  to  a 
material  point  subjected  to  strains  and  stresses  of  variable  amplitudes  which  produces  cracks  leading  to  total 
failure  after  a  certain  number  of  cycles"  and  it  is  the  main  cause  of  failure  of  machine  parts  in  service,  in 
mechanisms  and  structural  elements  in  aeronautics,  naval  and  automotive  industries.  Fatigue  failure  can 
occurs  under  load  conditions  well  below  the  strength  limit  of  the  material.  Typically,  a  progressive  loss  of 
strength  occurs  depending  on  the  number  of  stress/strain  cycles,  reversion  index,  load  amplitude,  etc.  This  loss 
of  strength  induce  the  material  to  inelastic  behaviour,  which  may  be  interpreted  as  micro-cracking  followed  by 
crack  coalescence  leading  to  the  final  collapse  of  structural  parts. 
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ORGANIZATION 


This  study  is  based  on  the  hypothesis  that  fatigue  damage  is  essentially  of  the  same  nature  as  mechanical 
damage  and  can  be  described  via  an  internal  variable  allowing  the  adequate  treatment  of  the  accumulation  and 
localization  of  dislocations,  therefore  the  theoretical  structure  of  continuum  mechanics  (such  as  plasticity  and 
damage)  is  suitable  for  the  study  of  non  linear  fatigue  problems.  The  inelastic  theories  of  plasticity  and/or 
damage  solve  the  problem  of  material  behaviour  beyond  the  elastic  range  and  both  theories  allow  the  study  of 
the  change  in  strength  that  a  material  point  suffers  by  inelastic  effects,  however  they  are  not  sensitive  to  cyclic 
load  effects.  In  this  work  the  standard  inelastic  theories  are  modified  to  account  for  fatigue  effect  coupled  with 
non-fatigue  behaviour  [1]. 

Concerning  composite  materials,  in  this  work,  a  numerical  model  to  assess  the  constitutive  non-linear 
behaviour  of  a  unidirectional  lamina  is  conceived  [2].  It  is  based  on  the  appropriate  management  of  the 
constitutive  models  of  component  phases  within  a  continuum  framework.  With  such  aim,  a  strategy  is 
developed  to  decouple  phases  that  may  also  consider  possible  inelastic  strains  and/or  elastic  degradation  of  the 
stiffness  based  on  previous  development  of  the  authors  [3],  [4].  Input  data  are  the  mechanical  properties  of 
each  component  (i.e.:  fibres  and  matrix  in  case  of  fibre-reinforced  laminates)  and  its  distribution  inside  the 
composite.  This  methodology  of  mixing  components  with  the  strategy  for  phases  decoupling  allow  each 
component  to  retain,  unchanged,  its  original  constitutive  law  (isotropic  or  anisotropic,  linear  or  non-linear) 
while  conditioning  the  composite  global  response.  Different  constitutive  models  (plasticity,  damage,  fatigue, 
creep,  etc.)  for  each  phase  may  be  taken  into  account. 

The  composite  model  (for  a  single  lamina)  is  combined  with  classical  lamination  theory  to  describe  laminates 
consisting  of  unidirectional  continuously  reinforced  layers.  Its  relative  simplicity  and  the  resulting  numerical 
efficiency  make  this  approach  well  suited  for  implementation  as  a  material  model  in  Finite  Element  programs 
for  studying  the  elastoplastic  response  of  structures  or  components  made  of  continuously  reinforced  or 
laminated  composites.  The  previous  mentioned  fatigue  model  is  applied  to  each  component  in  order  to  obtain 
the  durability  of  the  whole  composite  laminate. 


2.0  CONTINUUM  MECHANICAL  MODEL  FOR  FATIGUE  ANALYSIS 

It  is  assumed  that  each  point  of  the  solid  follows  a  damage-elasto-plastic  constitutive  law  with  the  stress  (5) 
evolution  depending  on  the  free  elastic  strain  variable  {E")  and  a  set  of  internal  plastic  and  damage  variables 
q={d’,d}  =  {E'^,  E,d=i^  },  were  E’’  and  represent  the  plastic  part  of  the  strain  and  a  unit 

normalized  dissipation  composed  by  the  plastic  plus  damage  parts,  respectively. 

The  free  energy  for  isothermal,  isentropic  and  adiabatic  processes  and  for  small  elastic  strains  and  large  plastic 
strains  is  written  in  the  reference  configuration,  accepting  the  additivity  of  its  elastic  and  plastic  parts, 
as 


T  (E,^ ,  r/)  +  T  ^  (a  f  )=(1  -  r/ )  /  C  T  ^  (a  f )  (1) 

2m 

where  E  is  elastic  Green  strain  tensor,  is  the  material  density,  d  =  E  is  the  internal  mechanical  damage 
variable  for  the  damage  processes  d'"'  <  (^/  =  1  with  the  initial  value  d’"‘  =  k“'  provided  by  the  defect 

level,  and  the  initial  constitutive  tensor.  The  stress  tensor  in  the  reference  configuration  S.j  can  be 
expressed  as 
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(2) 

dEij 

For  the  plastic  behaviour,  the  general  forms  of  the  yield  F  and  potential  G  plastic  functions  take  into  account 
the  influence  of  the  current  stress  state,  the  internal  plastic  variables,  and  other  variables  such  as  the  number  of 
cycles  Fl\ 


)  =  g{Sy)  =  constant. 

where  f{Sg)  and  g{Sy)  are  the  uniaxial  equivalent  stress  functions,  K(Sy,/c,N)is  the  strength  threshold 
(see  Figure  1).  All  the  internal  variables  at  current  time  t  are  obtained  by  means  of  an  integration  process 
af  =  [  dfdt ,  starting  from  its  evolution  law  af  -iuMf  ,  where  A.  is  the  plastic  consistent  factor. 

JO 


The  damage  function  is  defined  as 


G^{Sy,K)  =  S{Sy)-K{Sy,K,N)  (4) 

where  S(Sy)  is  the  uniaxial  equivalent  stress  function  in  the  undamaged  space,  K{Sy,K,N)  is  the  same 
strength  threshold  as  in  (3)  and  k‘'  -d  -^^d  dt  the  damage  internal  variable  with  an  evolution  defined  as 
d  =  where  p  is  the  consistency  damage  factor.  In  both,  (3)  and  (4),  the  normalized  dissipation 

is  defined  as  k  =  +  k'^  =  )/ ,  where  , 2™“  are  the  Clausius-Duhem  dissipation  for  the 

current  plastic,  damage  process  and  its  maximum  capacity  of  the  solid  dissipation  at  each  point,  respectively. 

The  effect  of  the  number  of  cycles  on  the  plastic  and/or  damage  consistency  conditions  ( A  =  0 ,  G^=  0 )  is 
introduced  as  follows. 
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f{S,j)-K{S,j,K).f^^,{N,S^^„R)^0 

K{Sy,R,N) 

(5) 

S{Sy)-K{Sy,K).f^^AN,S^,„R)^0 

K{Sy\R,N) 

(6) 

where  0  <  <  1  represent  the  unit  normalized  reduetion  part  of  the  strength  threshold  K  — plastie  and/or 

damage  strength  evolution—  by  load  eyclic  effeet. 

S-N  Curves 

Stress-N°  of  cyeles  {S-N)  eurves  are  experimentally  obtained  by  subjecting  identical  smooth  specimens  to 
cyclic  harmonic  stresses  and  establishing  their  life  span  measured  in  number  of  cycles.  The  curves  depend  on 
the  level  of  the  maximum  applied  stress  and  on  the  ratio  between  the  lowest  and  the  highest  stresses 
{R-S^,„  ! In  previous  works  ([1],[5],[6])  an  exponential  function  to  approach  steel  and  aluminium 
experimental  S-N  curves  was  proposed.  This  function  depends  on  and  is  capable  of  dealing  with  any  value  of 
the  ratio  between  minimum  and  maximum  stress.  However,  it  is  a  bit  difficult  to  adjust  its  parameters  to 
obtain  a  good  approximation  to  experimental  curves,  which  are  usually  not  defined  over  the  whole  life-span  of 
the  material.  Instead  of  using  such  an  exponential  function,  here  the  experimental  data  is  introduced  by  points 
in  a  table-like  mode. 

Usually,  S-N  curves  are  obtained  for  a  fully  reversed  stress  state  {R-  S^^^  / =  -1 )  by  rotating  bending 
fatigue  tests.  This  zero  mean  stress  is  however  not  typical  of  real  industrial  components  working  under  cyclic 
loads.  Based  on  the  actual  value  of  the  R  ratio  and  a  basic  value  of  the  endurance  stress  Se  (for  R  =  -1)  the 
model  proposed  postulates  a  threshold  stress  5^^ .  The  meaning  of  is  that  of  an  endurance  stress  limit  for 

a  given  value  of  /?  =  S^^^  /  S^^ ;  if  the  actual  value  of  is  /?=  -1  then,  . 

S,,  =  -  5,)  •  (0.5  +  0.5  •  - > abs{R)  <  1 

=S^+{S^-S^)-  (0.5  +  0.5  / - >abs{R)  >  1 

STHRl  and  STHR2  are  material  parameters  that  need  to  be  adjusted  according  to  experimental  tests. 

Cyclic  Strength  Reduction  Function 

The  S-N  curves  proposed  in  the  previous  section  are  fatigue  life  estimators  for  a  material  point  with  a  fixed 
maximum  stress  and  a  given  ratio  R.  If,  after  a  number  of  cycles  lower  than  the  cycles  to  failure,  the  constant 
amplitude  cyclic  loads  giving  that  maximum  stress  (and  ratio  R)  are  removed,  some  change  in  is 

expected  due  to  accumulation  of  fatigue  cycles.  In  order  to  describe  the  variation  of  the  following 
function  is  proposed: 


fred{R,Ncycles)  =  exp(-50  •  (loglO(A/cyc/ex))^™^) 
BO  =  -logl0(5_  /5J/(logl0(/V^))"^^"^ 
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BETAF  is  a  material  parameter  and  Np.  the  number  of  eyeles  to  failure. 


Time  Advancing  Strategy 

An  advantage  of  the  methodology  presented  consists  in  the  way  the  loading  is  applied.  In  a  mechanical 
problem  each  load  is  applied  in  two  intervals,  in  the  following  order  (see  Figure  2), 


-Tracing  load,  (described  by  “ai”  periods  on  Figure  2).  It  is  used  to  obtain  the  stress  ratio  R  -  5'^;^  ^  at 
each  integration  point,  following  the  load  path  during  several  cycles  until  the  R  relationship  tends  to  a 
constant  value.  This  occurs  when  the  following  norm  is  satisfied. 


GP 


pi+1  _  pi 
^GP  ^GP 


R 


i+l 

GP 


0 


(9) 


where  R^p  =  |  is  computed  at  each  Gauss  interpolation  point  for  the  load  increment  “i”. 


Figure  2:  Schematic  time  advancing  ioads  representation,  ai  describes  the  tracking  ioad 
domain  on  the  i  intervai.  bi  describes  the  enveioping  ioad  domain  on  the  same  i  intervai. 

-Enveloping  load,  (described  by  “bi”  periods  on  Figure  2).  After  the  first  tracing  load  interval  (ai),  the 
number  of  cycles  N  is  increased  while  keeping  constant  the  maximum  applied  load  (thick  line  in  Figure  2)  and 
the  stress  ratio  R.  In  this  new  load  interval,  the  variable  is  not  the  load  level  (kept  constant)  but  the  number  of 
cycles. 

This  two-stages  strategy  allows  a  very  fast  advance  in  the  time  loading.  A  new  interval  with  the  two  stages 
explained  should  be  added  for  each  change  in  the  loading  level. 
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3.0  CONSTITUTIVE  MODEL  FOR  FIBRE  REINFORCED  LAMINA 

In  composite  materials  each  component  makes  the  composite  behaviour  dependent  on  its  own  constitutive  law 
according  to  its  volumetric  participation.  In  addition,  the  components  morphological  distribution  inside  the 
composite  becomes  essential,  giving  better  properties  to  the  composite  material  than  its  component  parts. 
Classical  Mixing  Theory  (CMT),  also  known  as  Rule  of  Mixtures  (ROM),  takes  into  account  the  volumetric 
contribution  of  components  but  not  its  morphological  distribution.  Therefore,  CMT  can  be  only  useful  to 
obtain  basic  properties  of  composite  materials  assuming  parallel  behaviour  (components  with  same  strains  in 
all  directions).  This  hypothesis  is  a  strong  limitation  for  the  use  of  CMT  to  predict  the  behaviour  of  most 
composites  and  consequently  modifications  to  this  theory  are  always  required. 

In  the  late  1970s,  a  simple  micromechanically-based  constitutive  relation  for  composite  materials  has  been 
developed  [8].  This  model  consists  of  aligned  continuous  fibres  embedded  in  an  elastoplastic  matrix. 
Approximate  stress  and  strain  fields  are  prescribed  for  matrix  and  fibres,  strain  coupling  in  the  axial  direction 
and  stress  coupling  for  all  other  stress  and  strain  components  [9].  The  model  account  for  load  sharing  (fibres- 
matrix)  in  the  axial  direction,  but  fibres  do  not  interact  with  the  deformation  of  the  matrix  in  any  other 
direction,  leading  to  a  response  that  underestimate  the  transverse  and  shear  stiffness. 

To  overcome  these  limitations,  a  general  numerical  procedure  for  non-linear  constitutive  behavior  simulation 
of  the  unidirectional  laminas  is  developed.  This  model  is  able  to  capture  the  simultaneous  serial  and  parallel 
behaviours,  which  takes  place  inside  the  composite,  independently  of  loading  directions.  It  is  based  on  the 
appropriate  management  of  the  constitutive  models  of  component  phases  within  a  continuum  framework. 
With  such  an  aim,  a  strategy  is  developed  to  decouple  phases  that  may  also  contribute  possible  inelastic 
strains  and/or  elastic  degradation  of  the  stiffness.  Input  data  are  the  mechanical  properties  of  each  component 
(i.e.:  fibres  and  matrix  in  case  of  fibre-reinforced  laminates)  and  its  distribution  inside  the  composite.  This 
methodology  of  mixing  components  with  the  strategy  for  phase  decoupling  allows  each  component  to  retain, 
unchanged,  its  original  constitutive  law  (isotropic  or  anisotropic,  linear  or  non-linear)  while  conditioning  the 
composite  global  response.  Different  constitutive  models  (plasticity,  damage,  fatigue,  creep,  etc.)  for  each 
phase  may  be  taken  into  account. 

3,1  Definition  of  Serial  and  Parallel  Parts  of  Strain/Stress  Tensors 

The  SP  model  considers  that  in  a  certain  direction  (or  directions)  the  compounding  materials  behave  “in 
parallel”,  while  their  behaviour  is  “serial”  in  the  remaining  directions.  For  this  reason,  it  is  necessary  to  define 
and  split  the  serial  and  parallel  parts  of  strain  and  stress  tensors. 

Different  types  of  serial-parallel  mixtures  may  be  considered:  from  the  pure  parallel  (same  strain  field  for  all 
phases)  or  the  pure  serial  (same  stress  field  for  all  phases),  up  to  the  cases  of  parallel  behaviour  in  one 
direction  (fibre  reinforced  matrices)  or  in  a  plane  (biphasic  laminates).  This  work  is  focus  on  the 
unidirectional  fibres  case. 

Unidirectional  Fibres  Case 

Being  ,  the  director  versor  that  determines  the  parallel  behaviour  (fibre  direction),  and  being  and  the 
director  tensors,  orthogonal  to  ,  that  determines  the  directions  where  serial  behaviour  takes  place,  the 
parallel  projector  tensor  N^  may  be  defined  as  follows: 

Np=ej0ej  (10) 
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It  is  a  second  order  tensor  that  gives  the  parallel  (to  fibres)  projection  of  a  generic  vector  v : 

=  Np  •  V 

It  is  now  defined  the  4*-ordered  parallel  projector  tensor  Pp : 

Pp=Np0Np 

The  serial  projector  tensor  P^  is  evaluated  as  its  complement: 

P  =I-P 

I  g  I  Ip 

Both  tensors  allow  finding  the  parallel  part  of  the  strain  tensor,  Ep : 

£-p  =  Pp  :  £■ 


and  its  serial  part,  : 


By  this  mean,  the  strain  state  is  split  into  its  parallel  and  serial  parts: 

£  —  Ep  +  f  j 


The  stress  state  may  be  split  analogously: 


(7  —  <7  p  +  <7^ 


where: 


f7p  =  Pp  :  (T 
a  s  s'.  ^ 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


3,2  Hypothesis  for  the  Numerical  Modelling 

The  present  numerical  model  formulation  is  based  on  the  following  hypothesis: 

1)  Component  materials  have  the  same  strain  in  parallel  (fibre)  direction, 

2)  and  the  same  stress  in  serial  directions. 

3)  Composite  material  response  is  in  direct  relation  with  the  volume  fractions  of  compounding  materials. 
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4)  Homogeneous  distribution  of  phases  is  eonsidered. 

5)  Perfeet  bonding  between  eomponents  is  also  eonsidered. 


3,3  Equations  that  Govern  the  Serial-Parallel  Behaviour 


Constitutive  Equations  of  Compounding  Materials 

A  general  ease  is  considered,  with  materials  that  may  undergo  degradation  not  only  because  of  plastic  strains 
growth  but  also  because  of  deterioration  or  elastic  damage  in  their  stiffness: 


w  _  mr(m  ^  c"  P  m  f~\d\ 

(7  -  j  y  s,  s  5  ^  j 


(20) 

(21) 


As  an  example,  the  particular  case  of  additive  plasticity  is  shown: 

"'C 


Constitutive  equations  of  compounding  materials  may  be  rewritten  taking  into  account  the  serial  and  parallel 
split  of  strain  and  stress  tensors  -eqs.  (16)  and  (17)-: 


1 _ 

1 

1 - 

CL, 

b 

1 _ 

C^d 

^PP 

1 

u 

Lj 

bp 

_ 

C^d 

^SP 

c^d 
^SS  _ 

- 

c  P 

c^d 

^PP 

=  Pp 

:  "C"'  : 

P. 

c^d 

=  Pp 

;  : 

Pc 

where:  ■ 

^PS 

s 

with; 

c  = 

c^d 

^SP 

=  P. 

:  : 

Pp 

f - 

o 

=  P. 

:  : 

p. 

Equilibrium  and  Compatibility  Equations 

The  equations  that  define  the  stress  equilibrium  and  establish  the  strain  compatibility  between  components 
arise  from  the  analysis  of  the  previously  formulated  hypotheses: 

Parallel  behaviour: 


£p 


£  P  —  £  p 


<Jp 


+  <jp 


(22) 

(23) 
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Serial  behaviour: 


_  mi  m  ^ 

S  ^  k  £  s 

where  "'k  and  ^  k  are  volume  fraetion  eoefficients,  fulfilling:  '”k  +  ^  k  =  \ . 

3,4  Exposition  of  the  Problem 

Given  the  composite  strain  state  s  at  the  present  time  (t  +  At),  and  known  the  components  internal 
variables  of  the  previous  step  (t),  find  the  strain  state  of  the  components  and  also  the  corresponding  stress 
states  at  the  present  time  that  fulfil  the  equilibrium,  compatibility  and  constitutive  equations  described  in  the 
previous  section. 

The  following  chart  point  up  the  known  and  unknown  variables  of  the  problem: 


= 


+  ^  k  E 


(24) 

(25) 


Known  variables: 

t+At 

w. 

free  variable 

' 

'[V]. 

internal  variables 

Unknown  variables: 

t+ht 

[-£]. 

dependent  variables 

t+M  r 

t+At 

’ 

[^-].  “'w. 

t+M 

updated  internal  variables 

m  / 

Variables  ^  and  ^  group  all  the  set  of  internal  variables  corresponding  to  the  components,  like  for 
example  internal  variables  of  damage  and/or  plasticity  that  define  the  materials  state: 


"cq 


3,5  Algorithm  for  Compounding  Behaviours 

It  is  now  introduced  the  algorithm  for  the  constitutive  model  of  a  composite  material  made  up  of  behaviours 
of  different  materials  (e.g.:  fibre  and  matrix).  It  will  be  seen  that  the  model  that  sets  out  constitutes,  in  fact,  a 
manager  of  non-linear  constitutive  models,  which  combines  the  hypotheses  enunciated  for  composites  with 
general  material  models  (simple  or  complex  ones). 

Each  increase  of  load  may  cause  evolution  of  the  internal  variables  of  the  component  materials,  producing  an 
stress  disequilibrium  incompatible  with  the  second  hypothesis  (see  section  3.2).  This  stress  disequilibrium 
constitutes  the  residue  to  minimize: 
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Ao-5  =  "(7s  -  ^ (7s  (26) 

It  is  chosen  as  unknown  the  total  serial  strain  of  one  of  the  materials  -in  this  case  the  one  from  the  matrix 
( Ss  )-  since  the  strain  of  the  other  material  depends  on  the  first  one: 


^S\  )  fk  fk 

A(75(“f5)  ="as('”Ss)  -^(7s[^ £s["‘£s)) 

The  derivative  of  the  objective  function  (residue)  with  respect  to  the  unknown: 

-  ^0-5]  ^  a  "'(7s  _  d  .  a  ^£s 
a "£s  a "£s  a a 'a 

gives  the  expression  of  the  Jacobian: 


T 

k 

d"£s 

d"£s 

d  ^£s 

•  S 

= 

'_"cy 

1 

1 - 1 

0 

1 _ 1 

(-iih 

- 

;"cl; 

‘ 

that  it  is  used  to  update  the  unknown  at  each  local  iterative  step  of  the  Newton-Raphson  method: 

[".,r =[■.,]* -j-':[A<x,r 

Note  that,  from  now  on,  when  [•]  is  shown,  it  should  be  understood  [•]  . 


(27) 

(28) 


(29) 


(30) 


(31) 


Step  1.  Initial  Approximation. 

The  initial  approximation  of  the  unknown  can  be  established  when  considering  that  the  strain  increment 
maintains  the  tangent  evolution  of  the  previous  step: 

•  First,  the  constitutive  models  of  each  material  are  evaluated  to  determine  their  constitutive  tangent 
tensors  corresponding  to  the  previous  step. 


f  “1 

t 

r  1 

L 

9 

,L 

9 

L  “J 

Matrix 

Fibre 

Model 

Model 
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•  Then,  the  inerement  of  the  unknown  is  determined  supposing  that  the  total  strain  inerement  is 
distributed  among  eonstituent  materials  according  to  their  previous  step  tangent  stiffness: 

[Affj]  =  “'[e,]  -  '[«,] 

[As,]  =  “'[s,]  -  ■["£,] 

[-As,]”=A:['Cj,:[As,]  +  '-«:(^C-  "C)  :[As,]' 

where:  A  =  (  "'A: 

— I  ^ 

“  Sg  is  established. 

["s,]‘='[-s,]  +["As,]" 

k  =  a 

Step  2.  Residue  Evaluation, 

The  total  strain  tensors  of  the  components  |  j  are  determined  as  a  function  of  the  updated 

value  of  j  . 


[-s]‘  =["s,]  +  ["s,]‘ 

where 

=  -[s,] 

where 

The  constitutive  models  of  each  material  are  checked  to  determine  their  updated  internal  variables  and  stress 
states. 


["'f 

'[ 

s 

'[ 

^  a 

<=  < 

<=  < 

Matrix 

Fibre 

Model 

Model 

u 

u 

T;  [' 

T ;  [' 

a] 

Thereafter,  the  residue  is  evaluated  to  find  out  if  the  stress  equilibrium  (i.e.,  the  convergence  of  the  model)  has 
been  achieved. 
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[A<x,r  = 


where: 


Step  3.  Convergence  Checking. 

In  order  to  choose  the  tolerance,  the  order  of  magnitude  of  the  serial  stresses  must  be  considered  as  a 
reference.  If  the  stresses  of  the  previous  step  are  different  to  zero,  it  is  taken  the  minimum  between  them;  on 
the  contrary,  the  linearized  stresses  are  taken  as  a  reference. 


re/1  =  min  III  ‘["(7  s] 
■e/2  =  mm{  ■["C;,]:[s,] 


1 

'["CL 

]4^4  } 

If  ref\  >  0  then;  refer  =  ref\ 
else;  refer  =  ref  2 


The  tolerance  is  chosen  based  on  the  reference  value: 


toler  =  refer  -10"^ 


If  the  norm  of  the  residue  is  greater  than  a  tolerance,  then  to  go  to  step  4  -correction  of  the  unknown-, 
otherwise  go  to  step  5  -update  of  variables-. 


If 


>  toler 


then:  goto  Step  4. 


else:  goto  Step  5. 


Step  4,  Correction  of  the  Unknown, 

Evaluation  of  tangent  constitutive  tensors  of  both  materials. 

The  constitutive  models  of  each  material  are  evaluated  to  determine  their  tangent  constitutive  tensors. 
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Calculation  of  the  Jacobian. 


II 

1 - 1 

n 

1 _ 1 

*+^[-^CF]*  where:  [ 

m^T  _  p  . 

^ss]  • 

["C^T  ^ 

[ 

'CF]‘=Pr: 

['Cf  :  P, 

Update  of  the  unknown. 

[".,]*  :=[■.,]* -J 
k  \=k  +  \ 

:  [Aer,]* 

Go  to  step  2  -residue  evaluation-. 

Step  5.  Update  of  Variables. 

Once  obtained  the  convergence  of  the  constitutive  model,  all  the  variables  of  the  model  must  be  updated: 

f+A/ r  m  ~\  V  m  m  ~\  V  m  t+iS.tV  m  _ ~\  F  m  _ ”1^ 

L  <^sJ  =  L  ;  [  aj  =  [  aj  ;  [  C7j  =  [  crj 


Step  6.  Update  of  the  Composite  Stress  State. 

The  stress  tensor  of  the  composite  material  is  calculated  taking  into  account  equations  (23),  (24)  and  (17): 

The  presented  algorithm  manages  to  compound  of  behaviours  of  materials  fulfilling  the  hypotheses  of  the 
series-parallel  problem.  Note  that  the  same  algorithm  may  also  deal  with  the  extreme  cases  of  pure  serial  or 
pure  parallel  behaviours. 

3.6  Tangent  Constitutive  Tensor  for  the  Composite 

Once  obtained  the  local  convergence  of  the  composite  model  (manager  of  constitutive  models),  it  is  necessary 
to  find  out  the  composite  tangent  constitutive  tensor  that  will  allow  to  achieve  the  convergence  of  the  global 
problem. 

In  order  to  get  this  tangent  tensor  it  is  necessary  to  differentiate  the  composite  stress  tensor  resulted  form  the 
solution  of  the  linearized  system  of  governing  equations. 

Note  that  in  the  linearization  it  is  used  the  tangent  constitutive  tensors  of  the  component  phases  “C^  and 
(provided  by  each  simple  model)  which  are  based  on  the  updated  strain  states  and  internal  variables 
corresponding  to  each  component. 


m 


/  t+At  r  m  ~\ 

t-kAt  r  m  ~]\ 

1  L  ^J’ 

L  «J] 

(32) 
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C.  j- 

=  g 


(33) 


The  expression  of  the  algorithmic  tangent  tensor  is  obtained  differentiating  the  composite  stress  tensor  with 
respect  to  strain  tensor  in  a  way  consistent  with  the  integration  algorithm: 


where: 


dA(jp 

5Acr^ 

d<j  dAa 

dASp 
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ds  dAs 

dAcTp 

5Acr^ 

dAsg 

dAs^ 

y-^PP 
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,  fr  f^T  \ 

pp  '  A,  pp  1 

+  '"k  ^k('"Cls 

^PS 

"’C^  +  ^k 

^SP 

=  ("'k 

+  ^k 

^SP  ^ 

.  k.  f 

'-'SP  ) 

^ss 

\  ^SS  ^ss 

with:  A  = 

(  "'k 

ss  ^  ^ 

(34) 


It  can  be  clearly  observed  that  the  expression  of  the  composite  tangent  tensor  maintains  its  symmetry 
when  the  tangent  tensors  of  components  ( and  )  are  also  symmetrical. 

3,7  Enrichment  of  the  Transversal  Behaviour 

The  assumption  of  equal  stress  in  orthogonal  directions  to  the  fibre  (pure  serial  behaviour  in  transverse 
directions),  constitute  a  lower  bound  for  the  transverse  stiffness  of  the  composite.  Experimental  results 
confirm  this  fact  (see  next  figure).  For  this  reason,  an  enrichment  in  the  proposed  numerical  model  is  required 
to  be  able  to  predict  the  transversal  behaviour  more  accurately.  That  is,  the  Basic  SP  model  is  upgraded  into 
the  Enriched  Serial-Parallel  model  (ESP  model). 

In  order  not  to  resort  to  a  more  complex  model,  an  improvement  is  performed  that  maintains  the  structure  of 
the  previously  developed  algorithm:  the  behaviour  of  the  most  deformable  material  (matrix)  is  stiffened  by 
means  of  a  gamma  parameter.  The  structure  of  the  new  algorithm  is  similar  to  the  previous  one,  being  its  only 

difference  the  use  of  instead  of  ""f  ,  “cr,  whose  expressions  are: 


(35) 

('”o-)*  =  '"K:'”o-  (36) 

=  "K:  :  “K  (37) 
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where:  "’K  =  +  '> 

By  means  of  miero-mechanieal  eonsiderations  it  is  possible  to  take  values  of  the  gamma  parameter  as  a 
funetion  of  the  fibre  volumetrie  fraetion  (  and  the  ratio  between  fibre  and  matrix  Young  modules 

{^R  =  ^ E !  "'E^  (see  Figure  3).  For  example:  =  1.31  for  i?  =  20  and  =  0.6 . 


Volumetric  fraction  of  fibers 


R=60 
-©— R=40 
— A— R=20 
— A— R=10 
— R=5 
— B— R=2.5 


Figure  3:  Estimative  vaiues  of  gamma  as  a  function  of  the  fibre 
voiumetric  fraction  for  different  stiffness  ratios. 

In  the  following,  the  predietive  eapaeity  of  pure  Serial-Parallel  (SP)  model  and  enriehed  SP  model  for  the 
transversal  stiffness  of  a  glass-epoxy  lamina  with  Ef/Em=21.19,  Vf=0.22,  Vm=0.38  (stiffness  ratio  and  Poisson 
ratios)  is  compared  against  experimental  data  and  against  the  estimation  given  by  broadly  used  semi-empirical 
formulas.  An  hexahedral  composite  element  is  subjected  to  pure  transversal  load,  at  different  fibre  volume 
fractions.  In  Figure  4,  the  adimensional  curves  E2/Em  vs.  resulting  from  numerical  simulations  are  given 
together  with  experimental  values  taken  from  Barbero  [10]  (pag.  72).  The  curves  resulting  from  perfect 
inverse  R.O.M.  and  from  Flalpin-Tsai  equation  are  also  given  in  the  same  figure. 
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ORGANIZATION 


Figure  4:  Relative  transversal  stiffness  of  a  glass-epoxy  lamina  as  a  function  of 
fibre  volume  fraction.  Numerical  results  compared  against  experimental  data 
and  against  estimations  given  by  broadly  used  semi-empirical  formulas. 


The  transversal  stiffness  obtained  by  pure  SP  model  results  to  be  slightly  greater  than  the  one  given  by  inverse 
ROM,  due  to  Poisson  effeets.  Of  eourse,  when  Poisson  ratios  are  set  equal  to  zero  a  perfeet  inverse  ROM 
eurve  is  reeovered. 

The  graph  shows  that  simple  SP  model,  as  well  as  inverse  ROM,  underestimates  the  experimental  values.  On 
the  other  side  enriehed  SP  model,  with  an  appropriate  gamma,  obtains  an  approximation  to  experimental  data 
as  good  as  the  one  given  by  Halpin-Tsai  equation. 


4.0  CALIBRATION  AND  VALIDATION  OF  THE  MODEL 

The  numerieal  model  proposed  to  simulate  the  behaviour  of  eomposite  materials  is  based  on  the  appropriate 
management  of  the  eonstitutive  models  of  the  different  eomponent  materials  eonsidering  the  morphology  and 
distribution  of  fibres  within  the  matrix  [2]. 

Some  meehanieal  parameters,  sueh  as  the  transverse  Young’s  modulus  for  anisotropie  fibres  (e.g.:  earbon 
fibres),  are  impossible  to  obtain  by  means  of  direet  measurements.  Therefore,  it  must  be  dedueed  indireetly  by 
means  of  ealibration  of  this  parameter  to  fit  the  stiffness  of  the  lamina,  or  by  means  of  a  formula  that  eonsiders 
the  transverse  Young  modulus  of  the  eomposite,  the  Young  modulus  of  the  matrix  and  the  volumetrie 
partieipation  of  the  fibre  [10].  Onee  these  parameters  are  determined  they  remain  linked  to  that  eomponent 
material  sinee  they  are  its  attribute;  they  are  independent  of  fibres  direetion,  their  volumetrie  partieipation,  the 
loads  or  the  kind  of  problem. 

In  order  to  study  the  fatigue  phenomenon  in  eomposed  materials,  it  is  neeessary  to  ealibrate  the  parameters 
that  determine  the  life  eurves  S-N  of  the  matrix  and  of  the  fibre,  with  the  tests  done  on  laminas  of  eomposite. 
Then,  the  tests  on  laminates  and  struetures  will  allow  verifying  if  the  ealibration  has  been  made  sueeessfully. 
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As  calibration  methodology  for  the  fatigue  phenomenon,  the  following  tips  set  out: 

•  calibrate  the  life  curve  S-N  of  the  matrix  with  fatigue  tests  done  on  laminas  with  fibres  at  90  degrees; 

•  calibrate  the  life  curve  S-N  of  the  fibre  with  fatigue  tests  done  on  laminas  with  fibres  at  0  degrees; 

•  verify  that  the  calibrations  respond  to  the  behaviour  of  laminas  at  90  degrees;  failing  that,  return  to  the 
first  step; 

•  verify  with  the  response  of  the  fatigue  tests  done  on  [-l-45/-45]sym  laminates;  failing  that,  it  will  be 
necessary  to  include  an  interphase  "material"  and  to  calibrate  it  (currently  the  model  does  not  include 
the  interface  since  perfect  adhesion  between  fibre  and  matrix  is  assumed,  but  this  one  could  be  taken 
into  account  by  means  of  slight  changes  in  the  strain  compatibility  equations); 

•  verify  the  values  with  other  laminates  and/or  structures. 

The  model  is  validated  with  experimental  testing  done  at  ‘Institute  for  Static  and  Dynamics’  (ISD)  of 
University  of  Stuttgart  on  CFRR  material  samples.  Carbon  fibre  here  is  Tenax  HTA  5241  and  the  resin  used 
is  Ruetapox  VE4434.  Values  for  mechanic  properties  of  components  have  been  given  by  manufacturers  and 
they  are  listed  in  Table  1. 


Table  1:  Materials  mechanical  properties  given  by  manufacturers 


Components  Properties 

Unit 

Fibres  values 

Resin  values 

Density  (me) 

kg/m3 

1770 

1092 

Tensile  strength 

MPa 

3950 

Tensile  strain 

% 

1.5 

Young’s  modulus  parallel 

GPa 

238 

4 

Young’s  modulus  perpendicular 

GPa 

28 

4 

Shear  modulus 

GPa 

50 

1.5 

Poisson’s  ratio  (Vxv) 

- 

0.23 

0.35 

Mass  (Me) 

g 

3.771 

2.381 

Volumetric  content  ratio  (kc) 

- 

0.4942 

0.5058 

The  lay-up  has  been  chosen  to  get  a  symmetric  laminate  configuration. 

•  3 -layer  specimen  [0  90  0]  [90  0  90]  [0  90  0] 

•  4-layer  specimen  [90  0  90]  [45  -45  45]  [45  -45  45]  [90  0  90] 

•  5-layer  specimen  [90  0  90]  [45  -45  45]  [90  0  90]  [45  -45  45]  [90  0  90] 

The  thickness  of  the  laminas  is  chosen  in  such  a  way  that  the  numbers  of  fibres  are  the  same  in  both 
directions,  inner  layers  twice  thick  as  the  outer  ones. 

Static  testing:  For  each  type  of  specimen  (3,  4,  5-layers)  five  force  controlled  tensile  test  (European  Standard 
EN  2561)  have  been  carried  out  to  find  out  the  Young’s  modulus,  the  Poisson’s  ratio  and  the  ultimate  strength 
(stress  at  load  to  fracture). 

Figures  5  and  6  show  the  comparison  between  experimental  results  and  numerical  simulations  for  the  3 -layer 
specimens  at  0°  and  45°  loading  directions.  Eooking  at  the  break  pattern  of  the  specimens  it  can  be  seen  that 
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failure  is  caused  by  fracture  of  fibres  in  the  0°  loading  direction  (quasi-linear  stress-strain  curve)  and  due  to 
shear  failure  of  the  matrix  in  the  45°  case  (strong  non-linear  stress-strain  curve).  Both  situations  are 
accurately  reproduced  by  the  numerical  simulations. 


TCI 


0  - ! - 1 - ! - 1 - 1 - 1 - 1 - 1 - 1 - 1 

0  0.001  0.002  0003  0004  0.005  0.006  0.CO7  0008  OOS  001 

Parallel  Strain 


Figure  5:  Tensile  Strength  Test  -  3-layer  specimen  [0  90  0]  [90  0  90]  [0  90  0]. 

Dynamic  testing:  Tension  cyclic  loading  with  a  ratio  Smin/Smax  =  0.1  and  maximum  stress  depending  on  the 
thickness  of  the  samples  have  been  used  for  durability  tests.  For  3-layered  specimens,  load  levels  of  440,  480 
and  500  MPa  are  applied  up  to  failure.  The  damage  progress  within  the  inspected  specimens  are  observed  in 
different  ways:  Stiffness  monitoring  (stiffness  decrease  is  a  concomitant  of  damage  progress),  hysteresis 
measurement  (dissipated  energy  equivalent  to  area  of  hysteresis  loop  and  damage  progress)  and  heat 
radiation.  In  Figure  7  (left  side)  resulting  stiffness  plots  over  the  number  of  load-cycles  for  different  load- 
levels  are  summarized.  On  the  right  side  of  Fig.  7  the  corresponding  results  of  numerical  simulations  are 
presented.  Good  agreement  is  found  for  the  440  MPa  load  case,  acceptable  for  the  480  MPa  case  and  it  is 
difficult  to  assert  about  the  500  MPa  loading  case.  Looking  at  the  experimental  results  for  the  three  different 
load-levels  a  significantly  increasing  variation  of  the  number  of  load-cycles  to  failure  can  be  observed. 
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Parallel  Strain 


Figure  6:  Tensile  Strength  Test  -  3-layer  specimen  [45  -45  45]  [-45  45  -45]  [45  -45  45], 


Figure  7:  Fatigue  behaviour  of  3-layer  specimens.  Experimental 
results  on  the  left,  numerical  simulations  on  the  right. 
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ORGANIZATION 


5.0  CONCLUSIONS 

The  Serial-Parallel  (SP)  model  is  combined  with  classical  lamination  theory  to  describe  laminates  consisting 
of  unidirectional  continuously  reinforced  layers.  Its  relative  simplicity  and  the  resulting  numerical  efficiency 
make  the  SP  approach  well  suited  for  implementation  as  a  material  model  in  Finite  Element  programs  for 
studying  the  elastoplastic  response  of  structures  or  components  made  of  continuously  reinforced  or  laminated 
composites.  In  addition,  it  requires  relatively  small  computational  resources  when  implemented  into  a 
structural  FE  code.  Its  initial  drawback  of  underestimation  of  the  transverse  and  shear  stiffness,  is  then 
improved  upon  with  the  Enriched  SP  model. 

A  fatigue  model  based  on  continuum  mechanics  accounting  for  elasto-plastic-damage  constitutive  equations, 
previously  developed  by  the  authors,  has  been  applied  to  the  components  of  CFRR  material  in  order  to  obtain 
the  durability  of  the  whole  composite  laminate.  Comparison  between  experimental  and  numerical  testing 
carried  out  on  material  samples  allow  saying  the  methodology  presented  here  is  very  promising  for  durability 
analysis  of  composite  materials  and  structures. 
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SYMPOSIA  DISCUSSION  -  PAPER  NO:  31 

Author’s  name:  O.  Salomon 


Discussor’s  name:  J.  Calcaterra 

Question:  In  the  continuum  fatigue  formulation,  you  assumed  a  fatigue  endurance  limit  stress-  Are  there  data 
to  validate  this  assumption  for  your  material? 

Answer:  For  composite  materials  experimental  data  corresponding  to  components  must  be  obtained  from  tests 
on  the  final  material  (composite)  using  special  configurations  and  ‘known’  data  for  one  of  the  composites. 
This  is  not  an  easy  task  but  can  be  done.  Anyway,  a  ‘fixed’  stress  endurance  limit  is  not  a  requirement  for  the 
model.  If  stress  continues  decreasing  as  a  function  of  number  of  cycles,  this  information  is  the  important. 
There  is  no  need  for  a  constant  endurance  limit. 
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